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Abstract. The Cosmic Neutrino Background (C^B) anisotropy is calculated for 
massive neutrino states by solving the full Boltzmann equation. The effect of weak 
gravitational lensing, including the Limber approximation, is also derived for massive 
particles, and subsequently applied to the case of massive neutrinos. 



1. Introduction 

The Cosmic photon Microwave Background (CMB) is currently our main source of 
information about the physical content of the Universe. Observations of the CMB 
anisotropy provides detailed information about the curvature of the Universe, the matter 
content, and a plethora of other parameters [I]. 

Standard model physics likewise predicts the presence of a Cosmic Neutrino 
Background (Cz/B) with a well defined temperature of T„ ~ (4/H) 1 / 3 ^. While it 
remains undetected in direct experiments, the presence of the Cz/B is strongly hinted 
at in CMB data. The homogeneous Cz/B component has been detected at the 4-5er 
level in the WMAP data (see e.g. [TJ El EH HI El Ej). Furthermore, this component is 
known to be free-streaming, i.e. to have an anisotropic stress component consistent 
with what is expected from standard model neutrinos (see [3, El El QUI EEU EEJ fl3]). 
Finally the standard model neutrino decoupling history is also confirmed by Big Bang 
Nucleosynthesis (BBN), the outcome of which depends on both the energy density and 
flavour composition of the Cz/B. 

While this indirect evidence for the presence of a Cz/B is important, a direct 
detection remains an intriguing, but almost impossible goal. The most credible proposed 
method is to look for a peak in beta decay spectra related to neutrino absorption 
from the Cz/B [T4"l [T5l [16], although many other possibilities have been discussed 
[mimiiaEfflEHESESEU. The neutrino absorption method was first investigated 
by Weinberg [H], based on the possibility that the primordial neutrino density could 
be orders of magnitude higher than normally assumed due to the presence of a large 
chemical potential. Although a large chemical potential has been ruled out because 
it is in conflict with BBN and CMB [25l ESJ E3 EH1 EH], the method may still work 
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and recently there has been renewed interest in detecting the Cz/B using beta unstable 
nuclei. 

Although the direct detection of the Cz/B is already very challenging, one might 
speculate on the possibility that in the more distant future anisotropies in the Cz/B 
will be detectable. For massless neutrinos the calculation of the Cz/B proceeds in a 
way which is almost identical to the standard CMB calculation. The massless Cz/B 
anisotropy spectrum was first presented in [3U] and subsequently calculated in [3]] in a 
highly simplified way which contains some, but not all of the essential physics. 

For massive neutrinos the calculation is much more complicated, and the Cz/B 
anisotropy is changed considerably: If the mass is sufficiently high, neutrino velocities 
can be as low as the escape velocities of galaxies. In this case the Cz/B is entirely 
determined by non-linear gravitational clustering. The current thermal velocity 
of a non-relativistic homogeneous neutrino background is given roughly by (v) ~ 
1500 km s" 1 ^ ° m° V J > which should be compared to gravitational streaming velocities 
which are up to ~ 1000 km s" 1 . For small masses (i.e. non-degenerate, m v < 0.1 eV) 
it is possible to make a calculation which is analogous to what is done for the CMB. 
As will be explained later the Cz/B spectrum is shifted to larger angular scales, mainly 
because the much shorter distance to the neutrino last scattering surface changes the 
relation between angular scales and length scales. Furthermore the amplitude of the 
anisotropy greatly increases at small multipoles because the gravitational source term 
becomes much more important as neutrinos become increasingly non-relativistic at late 
times. 

In addition to this change in the primary Cz/B spectrum, the effect of gravitational 
lensing is also very different from the case of massless neutrinos. As is the case for the 
primary spectrum, gravitational lensing becomes increasingly important at low / as the 
mass increases. A detailed calculation of the lensing of massive neutrinos is presented 
in section 3. First, however, we derive the necessary equations for the primary Cz/B 
anisotropy in section 2 and present a numerical calculation of the CuB power spectra 
for different masses. In section 4 we combine the results of sections 2 and 3 to derive the 
form of the lensed massive CuB. Finally, section 5 contains a discussion of our results 
and our conclusions. 

2. The primary Cz/B 

2.1. Theory - The Boltzmann equation 

The evolution of any given particle species can be described via the Boltzmann equation. 
Our notation is similar to that of Ma and Bertschinger |32j . As the time variable we 
use conformal time, defined as dr\ = dt/a(t), where a(t) is the scale factor and t is 
cosmic time. Also, as the momentum variable we shall use the comoving momentum 
qj = apj. We further parametrise qj as qj = qnj, where q is the magnitude of the 
comoving momentum and nj is a unit 3-vector specifying direction. 
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The Boltzmann equation can generically be written as 

L[f] = ^ = C[f], (1) 

where L[f] is the Liouville operator. The collision operator C[f] on the right-hand side 
describes any possible collisional interactions. For neutrinos C[f] = after neutrino 
decoupling at T ~ 2 — 3 MeV. 

We then write the distribution function as 

f{x\q,n j ,T 1 ) = f (q)[l + V(x\ q,nj,r))}, (2) 

where fo(q) is the unperturbed distribution function. For a fermion which decouples 
while relativistic, this distribution function is 

f (q) = [exp(g/T ) + l]" 1 , (3) 

where T is the present-day temperature of the species. 

In conformal Newtonian (longitudinal) gauge the Boltzmann equation for neutrinos 
can be written as an evolution equation for \l/ in /c-space 



jo or e dmq 



4> — ik-\ii\) 



0, (4) 



where \i = nPkj. ip and <fi are the metric perturbations, defined from the perturbed 
space-time metric in the conformal Newtonian gauge [32] 

ds 2 = a 2 (r/)[-(l + 2ifj)dr] 2 + (1 - 2<f>)5 ij dx i dx i ]. (5) 

The perturbation to the distribution function can be expanded as follows 

oo 

f^H/i^l)^^). (6) 

1=0 

One can then write the collisionless Boltzmann equation as a moment hierarchy for the 
\EVs by performing the angular integration of L[f] 

*o= -k-^ +<!>—, (7) 

* 1 = k±& -2* 2 )-±^, (8) 
3e 3q dmq 

= fc ( 2 / + l)6 (/ ^" 1 " (/ + ' 1 ~ ^ (9) 

By integrating the neutrino perturbation over momentum 
j dqq 2 ef Q {q)^ l 
J dqq 2 efo{q) 

one finds a set of equations equivalent to those used to follow the perturbations in 
photons or massless neutrinos, i.e a set of equations for F v \. 

The distortion to the sky intensity can also be found using the perturbation to the 
temperature, 9, related to the distribution function via 

/(g) = [exp(g/[T o (l + 0)] + l]- 1 . (11) 
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Equating Eqs. (jzj) and (JTTj) . ^ and 6 are related by 

Since the transformation between \1/ and is mass-independent, the normalisation of 
the angular power spectrum is not significantly affected. [3T] used the transformation 
— (-3j ) which introduces an extra factor of v 2 relative to our definition. The v 2 
factor, by construction, suppresses the perturbations for the higher neutrino masses 
significantly. 

By substituting into the Boltzmann equation, Eq. (J3J), it can be seen that 
is g-independent in the massless case. For massless particles it is therefore convenient 
to calculate the angular power spectrum of the temperature perturbation. For massive 
particles is g-dependent giving rise to spectral distortions in the temperature field. 
Since it is convenient to have results similar to the CMB in the massless neutrino case 
we calculate the massive neutrino anisotropy spectrum of the quantity Qi(q) related to 
Qui by 

f dqqhfoiq^q) 

Ul ~ fdqq 2 ef (q) ' [U) 
Of course, the quantity which is actually measurable will depend on the type of 
experiment used. A typical experiment will measure either a number flux, a momentum 
flux, or a kinetic energy flux as a function of angle. For massless particles these 
quantities are trivially related by a momentum independent number, i.e. the number 
flux anisotropy is 30„/ and the momentum/energy flux anisotropy is F v i = 40„/. For 
massive particles this is no longer true and one must calculate the appropriate quantity 
for any given experiment. 

Similarly to the photon case one can then construct the CVB sky brightness 
fluctuation angular power spectrum as 

Cf(q) = (4vr) 2 J k 2 dkPj(k)Q 2 (q, k). (14) 

Here Pi(k) is the primordial potential fluctuation power spectrum, Pi(k) oc k n ~ A . 
Throughout the paper we assume a flat Harrison-Zel'dovich spectrum so that n — 1 
with cosmological parameters (f2 b , fi m , h, A s ) = (0.05, 0.3, 0.7, 0.7, 2.3 ■ 10~ 9 ). 
As can be seen, a given / gets contributions from all k. The extra g-dependence arises 
because in principle one should perform the g-dependent lensing of Cf(q) in g-bins. The 
total Cf is found by averaging over momenta at the present time 770 

/rfgg 2 6(r / o,g)/o(g)y / Cp(gT j^ 
f dqq 2 e(r] ,q)f (q) 

Cf does not include lensing, since the second-order term encoding deflections has been 
left out of the Boltzmann equation. Note that mapping the neutrino anisotropic sky 
in different momentum bins (anisotropic neutrino momentum tomography) will probe 
structures at different spacial distances from us. Furthermore, the observed CuB will 
be a superposition of the spectra for each individual neutrino mass. 



Cf 
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Figure 1. Primary Cz/B spectrum for different neutrino masses. 



2.2. Gauge effects 

As noted we use the conformal Newtonian gauge, since this gauge is directly related to 
physically measurable quantities. In the Synchronous gauge the velocity perturbation, 
6, for the CDM component is, by definition, zero. Therefore 9 V , which is a momentum 
integral over is a gauge dependent quantity. In contrast, the anisotropic stress, 
which is a momentum integral over \I/ 2 , is a gauge independent quantity. Since all 
moments with / > 2 is recursively related to ^2, these higher order moments are 
gauge independent as well. 

Since we cannot separate the CMB/CVB dipole from our own peculiar motion, we 
are only interested in modelling the C®'s with / > 2 when comparing with observations. 
But C*i, the lowest mode containing physically relevant information, is gauge dependent. 
We have taken this into account by working in the physical conformal Newtonian gauge. 

We also note that the transfer functions are gauge dependent, though for the 
Synchronous and conformal Newtonian gauges they are almost identical inside the 
horizon for the massive components. Therefore we have calculated the transfer functions 
used to get the lensing contribution in the Synchronous gauge with CAMB [33] . 

2.3. Numerical results 

We have used the COSMICS code [31] to solve the Boltzmann hierarchy for the 
neutrinos. In practise we have solved the system going up to / = 500 with 64 bins 
in q, equally spaced from q/T = to 15. In Fig. [I] we show results for Cf for various 
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Figure 2. Sky maps of the primary neutrino power spectra, Cf , with the dipolc 
included, for m v = 1(T 5 eV (top-left), 1CT 3 eV (top-right), 1CT 2 eV (bottom-left) and 
10 _1 eV (bottom-right). The maps have been generated with the same underlying 
random numbers with the HEALPIX package [35] . 



masses and Fig. [2] shows sky map realisations for these spectra. 

The massless case (i.e. 10 _5 eV) is consistent with the result of [30]. At high I 
the spectra are almost identical, and do not depend on the neutrino mass. The reason 
for this can be understood from the following argument: Above a certain k- value, fcps, 
neutrinos are completely dominated by free-streaming and this A;- value is proportional to 
m„. In order to convert this to an /-value one then uses the relation ips ~ ^fsX* (where 
X* is the comoving coordinate from which the neutrinos originate) and since x* m u 1 
for non-relativistic particles [SSj, hs does not depend on m v . Inserting numbers one 
finds Zfs ~ 100 which is in good agreement with Fig. [TJ At smaller angular scales, 
I > l-ps, the anisotropy comes from the Sachs- Wolfe effect during radiation domination. 

For smaller /-values the anisotropy increases dramatically as the mass increases. 
This can be understood as follows. As soon as neutrinos go non-relativistic the ^ip d d\Jq 
term in begins to dominate the Boltzmann hierarchy evolution. This quickly makes 
the higher / modes increase as well, and the final amplitude simply depends on the time 
elapsed after neutrinos go non-relativistic. 

The effect can be seen in Fig. [3] which shows the evolution of ^f\, ^2 and \l/io for 
three different neutrino masses and two different /c-values. As soon as neutrinos go 
non-relativistic immediately begins to grow, and the higher \l/;'s follow with a slight 
delay for k = 0.1 h Mpc^ 1 . This exactly matches the low / behaviour seen in Fig. [TJ 




Figure 3. ^'s for 3 neutrino masses with momentum q/To = 3 as a function of the 
scale factor. The upper three panels are for k = 0.01 /iMpc -1 and the lower three 
panels for k = 0.1/iMpc -1 . 
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Note also the behaviour of \l/io for k = 0.01 /iMpc" and m„ = 0.1 eV. For this mass 
and k- value, the q/e term in the equations becomes sufficiently important to suppress 
the propagation of the gravitational source term to high I. At higher k this is no longer 
true. 

For high mass neutrinos (> O.leV) gravitational distortion is so strong that the 
dominant contribution comes from the local galactic and/or cluster halo. This case 
is similar to the study of the WIMP flux anisotropy and requires iV-body simulations 
[37] . This will be treated separately in a later paper while in the present paper we limit 
ourselves to the framework of linear theory. 



3. The lensing distortion 

3.1. Theory 

For strictly massless neutrinos the lensing distortion is identical to that for photons 
[38l l39l HO]. The change in angle da per unit path d\ is proportional to the transverse 
derivative of the gravitational potential, ip, 

da , , 

T X = - 2V ^ < 16) 

where V±ip is the component perpendicular to the line of sight. Relaxing the assumption 
of relativistic particles and solving the geodesic equation for a neutrino propagating in 
a weak potential, one arrives at the result 

_ = -_V^. (17) 

Note that this result reduces to the ordinary Newtonian expression ^ = — ^Vj,^ in the 
limit v — > 0. It should also be noted that the expression diverges as v — > because the 
assumption v > v csc is violated, i.e. particles with low velocity will be gravitationally 
bound in the potential. 

From Eq. ( TT7|) we can calculate the distortion spectrum in a manner similar to what 
is done for the usual gravitational lensing spectrum. This is done with the assumption 
of Gaussian perturbations and using the Born approximation [ID]. 

The total deflection angle a is related to the lensing potential, II, by 

a = VnlT. (18) 

In a flat universe the normal formula for the angular lensing power spectrum, 
(Ui m Il*, ml ) = SwSmm'Cj 1 , is given by 

* 2 

c? = i ^jf £ dxv^ik,^ - xMkx) (^r) ' (19) 

derived using Eq. (Tl6|) . Here the power spectrum, V^, is related to the ordinary matter 
power spectrum in the density contrast by oc a P m I k . X* is the conformal distance 
at which the photons (or neutrinos) decoupled, taken to be a single source sphere, and 
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ji is a spherical Bessel function. ji(kx){x* ~ x)/{x*x) is an effective window function, 
which distributes power in /c-space along the particle trajectory to angular Z-space. 

However, several changes are necessary when particles are allowed to have mass. 
The relation d\ = —dr\ must be replaced with d\ = —vdr] (the minus sign accounts for 
the fact that time and space run in different directions, i.e. J* d\ ~ P Q dr\ with rf ~ 0, 
so that the observer is at the origin). In addition, the power spectrum V^(k,r] — x) 
should be replaced by V^(k,rj). With these modifications, the expression for massive 
particles becomes 



Q n (s)=4vr|f[A?(g, 



with 



and 



A?(g,*) 



'A) 



drjAY(q,k,rj), 



x*(q) -x(q,v) 
x*(q)x(q,v) . 



(20) 



(21) 



(22) 



where v(q,rj) = q/e = 1/y 1 + a 2 (rj)m 2 / q 2 . x* 1S now momentum dependent since 
neutrinos with different velocities have different distances to their respective last 
scattering surfaces, though they still scattered at the same time rj*. 

We calculate an average quantity of the lensing power spectrum found by doing an 
energy average over the C z n (q , )'s 



J^g 2 e(r /0 ,g)/ (g)v / CfM 



(23) 



/ dqq 2 e( y r] ,q)fo{q) 

Using the orthogonality of the Bessel functions together with the fact that they 
pick out the scale k ~ l/x at high /, the high / limit of the above equations reduce to 
the Limber approximation (see e.g. [5H S2]) for massive particles 

2tt 2 rvo 



Cf{q) 



P 



v 



x -x 
x*x 



(24) 



For I > 100 this approximation is very good for all masses simulated. 



3.2. Numerical results 

In Fig. H] we show/^Z + l) 2 ^ for various neutrino masses, all in a ACDM background 
model, and Fig. [5] shows realisations of these spectra. We have done the calculation in 
linear theory only. As can be seen, for higher neutrino masses the lensing distortion 
peaks at lower / because a given /c-scale corresponds to lower I when v < c. Basically 
there is no contribution from modes with k < l/x*- This also means that for higher 
masses Cj 1 picks up a much larger contribution from high k. This can be seen explicitly 

X Compared to Cf their is an extra factor of 1(1 + 1) since the physically relevant quantity is the 
deflection angle. 
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Figure 4. The lensing potential angular power spectrum, Cj 1 , for 5 different neutrino 
masses. 




Figure 5. Sky maps of the lensing deflection 1(1 + l)C; n with the dipole included, 
for m v = 1Q- 5 eV (top-left), 1CT 3 eV (top-right), 1CT 2 eV (bottom-left) and 1CT 1 eV 
(bottom- right). 
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Figure 6. The contribution to Cj 1 from various scales k in units of hMpc . From 
top-left to bottom-right the neutrino mass is 10~ 5 eV, 10 -3 eV, 10 -2 eV and 10 _1 eV. 

in Fig. [U] which shows the contribution to Cj 1 from various scales. For the higher masses 

there are significant high-fc contributions to lensing already at low I. 

From the Limber approximation it is also straightforward to understand how Cj 1 

changes with neutrino mass. At high /, very approximately we can set v oc \ lTy Eq. (124j) . 

r » i 2 

which means that ^ \*x 00 ^1^' ^' ie potential power spectrum changes from 
Vip(l/x, v) ~ Const, at low / to V^iljx-, v) ~ X 4 a ^ high /. Thus, at high / the integrand 
is proportional to \ 2 which is proportional to m~ 2 . This explains the lower overall 
lensing power at high / for high masses. 

From Figs. Q] and H] it can be seen that there is a large cross-correlation between 
Cf and Cj at low /. This was also shown in [40J. 

In Fig. [7] we show the contribution to lensing from different redshifts. Almost all of 
the low / contribution comes at very low z. Note that for the m v = 10~ 2 eV case which 
was semi-relativistic until a fairly low redshift the high-z contribution is fairly similar 
to the massless case, i.e. the transition from relativistic to non-relativistic can be seen 
directly from the change in shape of the lensing spectrum. 
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Figure 7. Contributions to the lensing potential, C { P, at different redshifts. From 
top-left to bottom-right the neutrino mass is 10 -5 eV, 10 -3 eV, 10 -2 eV and 10 _1 eV. 



4. The lensed CuB 



4-1. Theory 



Finally, in this section we combine the results from the two previous sections to derive 
the lensed CzaB spectrum. 

Modelling weak gravitational lensing as a second-order effect, the full-sky lensed 
angular neutrino power spectrum is found from 



Cf{q) = 2tt J ' i(q, cos/3, (25) 

where /?) is the lensed correlation function and cos/3 = fii ■ fi 2 , where fii and n 2 
indicate two directions on the sky. d l 00 ((3) is a special case of the reduced Wigner 
functions, d l mm ,((3), given by 

d l mm >(f3) = (-l) l - m '[(l + m)\{l - m)\(l + m')\(l - m')!] 1/2 (26) 
. k [cos(/3/2)] m+m ' +2fc [sin(/3/2)] 2Z - m " m '- 2fe 



k\(l — m — k)\(l — m! — k)\(m + m! + k)\ 
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Figure 8. Percentage difference between the iensed and uniensed angular power 
spectra, {Cf — C[ ■)/C' l ■ , for 2 different neutrino masses and compared with the lensing 
effect on the CMB. The spectra have been slightly smoothed due to finite numerical 
resolution. 



where the sum is over all k fulfilling the criterion that the arguments of the factorials 
should be non-negative. 

Taking sky curvature into account £(g, (3) is given by 

® - E ^^)e- L2c + M)/2 (27) 
i=i 

V d l mm ,({3) 1^ [L 2 C + (q, P)/2]I^ [L 2 C.(q, {3)/2}. 

L ' 2 2 

mm' 

Here Cf(q) is the uniensed power spectrum and in the double sum, J^ mm /, m and m! 
runs from — I to + I in integer steps with the criterion that m + m' is even. L = 1 + 1/2 
and /^(x) is the modified Bessel function of the first kind. Finally 

00 21 + 1 

c±(?, /*) = £ ^^(^ + l)C?(?)4u03). (28) 
(=1 

where C z n (g) is the power spectrum of the lensing potential. 
^.J?. Numerical results 

We have found C® with the lensing code in CAMB and Fig. [8] shows the difference with 
respect to the uniensed spectra for the m v = 10 -5 eV and 10~ 3 eV well as 

for the CMB. Both neutrino spectra are closely correlated with the CMB case because 
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the main effect is to move power from the higher I modes in the primary spectra which 
are very similar at these /-values. Since the primary spectrum dips to a minimum at a 
somewhat higher / for m v = 10~ 3 eV than for m v = 10~ 5 eV the lensing effect also kicks 
in later, as can be seen in Fig. [HJ 

We have not presented lensing results for the higher masses because the relative 
effect of lensing becomes less important at low I, and because numerical noise from the 
truncation of both primary and lensing spectra at I = 500 prevents a reliable calculation 
of the lensed CuB spectrum much beyond / ~ 10 for the most massive cases. 

Note that we have found Cf by lensing the average primary spectrum with the 
average lensing spectrum. For small neutrino masses with v ~ c this is surely a good 
approximation. For higher masses it will be more important to lense each Cf(q) with 
its corresponding Cj l (q) and then calculate an energy average, though the order of 
averaging should not significantly affect the total Cf. 

5. Discussion and conclusions 

We have calculated the anisotropy of the Cz/B in linear theory which applies to neutrino 
masses of less than ~ 0.1 eV. For massless neutrinos the power spectrum of CvB 
fluctuations closely resembles the usual CMB spectrum, but with the baryon-photon 
acoustic oscillations absent. 

At high / the neutrino spectra are almost identical, independent of the neutrino 
mass. The reason is that at high / all neutrinos are dominated by free-streaming which 
in /-space has approximately the same impact for all masses. 

For smaller /-values the anisotropy increases dramatically as the mass increases, 
because the gravitational source term becomes much more important at late times for 
massive particles. This initially increases the lowest multipoles but via the Boltzmann 
hierarchy the effect quickly propagates to higher /. 

We then proceeded to calculate the effect of weak gravitational lensing for massive 
neutrinos and found it to be much stronger at low /, and correspondingly weaker at high 
/, as compared to the massless case. Finally we calculated the effect of lensing on the 
primary CvB spectra and found the effect to be unimportant (with relative changes at 
the per mille level up to / ~ 50), but with some differences depending on the neutrino 
mass. 

It is worth mentioning that any direct experimental measurement of the Cz/B 
anisotropy will most likely measure flavour states, not mass states. The actual 
anisotropy measured will therefore be a superposition of anisotropics for three different 
mass states, weighed with their individual flavour content. 

We should finally again stress that our results are only valid for masses of < 0.1 eV. 
For higher masses linear perturbation theory breaks down because neutrino streaming 
velocities become comparable to the typical gravitational flow velocities so that a 
significant fraction of neutrinos are bound in structures. In this case the CvB spectrum 
must be found from iV-body simulations of neutrino clustering [37]. This can also be 
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seen from the fact that the anisotropy at low I is a factor ~ 10 9 higher for m u = 0.1 
eV than for massless neutrinos. Since the anisotropy for massless particles corresponds 
to 5p/p ~ 10~ 5 the corresponding Sp/p for 0.1 eV neutrinos is of order one, indicating 
that perturbation theory breaks down. 
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